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ABSTRACT 


This  report  describes  a  VHF/UHF  propagation  prediction  model  and  a  series 
of  field  measurements  aimed  at  verifying  the  results  of  the  model.  Comparisons 
between  the  model  and  the  trial  results  were  made  at  spot  frequencies  from  100 
MHz  to  900  MHz  for  several  types  of  antenna  and  for  various  heights  (not 
exceeding  2  metres)  above  ground  level. 

The  predicted  path  losses  from  the  model  were  found  to  agree  quite  well 
with  the  results  of  the  field  measurements.  It  should  be  noted  that  the 
prediction  model  was  never  Intended  to  be  a  comprehensive  tool,  merely  a 
utility  package  which  could  be  incorporated  into  other  software. 

The  program  was  written  on  an  HP  9000  series  model  310  computer  using 
BASIC  5.1.  This  language  has  a  set  of  commands  which  allow  direct  manipulation 
of  complex  values,  consequently,  translation  into  another  language  or  version 
of  BASIC  not  having  these  facilities  will  require  extra  programming. 
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SHORT  RANGE.  CLOSE  TO  GROUND.  VHF/UHF  PROPAGATION: 

A  COMPARISON  OF  MEASURED  RESULTS  WITH  A  SIMPLE  RAY-BASED  MODEL 


1.  Introduction. 

A  requirement  exists  to  study  the  propagation  of  signals  between  100  MHz 
and  1GHz  at  ranges  of  less  than  50m  and  with  antenna  heights  less  than  2m.  A 
propagation  loss  prediction  program  was  written  employing  geometric  ray- theory 
which  agreed  quite  closely  with  previously  published  data  (Ref.l).  However, 
this  data  did  not  include  results  obtained  when  one  or  both  of  the  antennas  is 
much  less  chan  lm  above  the  ground.  It  was  therefore  necessary  to  undertake  a 
aeries  of  measurements  using  various  types  of  antenna  to  validate  predictions 
obtained  from  the  model  under  the  conditions  for  which  no  data  was  available. 

This  report  describes  both  the  propagation  model  and  the  trials  devised  to 
validate  the  predicted  results . 


2.  Theoretical  background  to  the  prediction  model. 

Energy  radiated  from  a  transmitting  antenna  may  be  considered  to  consist 
of  three  parts  -  a  direct  wave  which  travels  along  a  direct  path  from 
transmitting  antenna  to  receiving  antenna,  an  earth- ref lected  wave  which 
arrives  at  the  receiving  antenna  following  a  reflection  at  the  earth’s  surface, 
and  a  surface  wave  which  Is  supported  by  and  travels  along  the  earth's  surface. 
The  direct  wave  and  the  ground  reflected  wave  (the  sum  of  which  is  known  as  the 
space  wave)  can  be  accounted  for  by  ray  theory  whereas  the  surface  wave 
requires  a  solution  of  suitable  electromagnetic  field  equations. 

The  propagation  model  calculates  propagation  loss  between  two  similar 
antennas  by  evaluating  the  sum  of  the  direct  ray  and  the  ground  reflected  ray 
at  the  receiving  antenna.  The  effect  of  the  surface  wave  Is  relatively  small  at 
frequencies  above  about  100  MHz  and  is  omitted  from  the  model  (it  does  become 
more  significant  if  the  antennas  heights  are  small  in  terms  of  wavelengths  and 
if  vertical  rather  chan  horizontal  polarisation  is  used  but  in  the  present 
context  the  error  is  still  small) . 


2.1  ground ..rofle<?U?n  sufficient. 

The  electrical  properties  of  the  ground  are  represented  by  the  following 
equation  for  refractive  index  (assuming  tr  -  1,  o  -  0,  for  air), 

n 


where  «r  -  ground  relative  permittivity 

a  -  ground  conductivity 
u  ™  2*f 

<Q  -  permittivity  of  free  space  (8.854X10  F/m) 
Mr  is  taken  as  unity  for  both  ground  and  air. 
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Typical  values  for  Che  frequency  range  under  consideration  are: 

For  wet  soil,  -  15  a  -  0.1 

For  dry  soil,  tr  -  4  a  -  0.0001 

The  following  equations  give  the  ground  reflection  coefficient  in  terms  of 
n2  (Ref. 3)  : 


For  vertical  polarisation, 

_  _  n2  sin  9  -  J( n2-cos2  0) 
n2  sin  t  +  V(n2-cos2  9) 


For  horizontal  polarisation, 


sin  9  -  y(n2-cos2  9) 
sin  9  +  7(n2-cos2  J) 


where  n2  is  defined  above  and  9  is  the  angle  of  incidence  between  the  ray 
and  the  horizontal . 

If  these  functions  are  represented  graphically  it  will  be  seen  that  for 
small  angles  of  incidence  the  reflection  coefficient  for  horizontal 
polarisation  is  approximately  -1.  With  vertical  polarisation  |R|  falls  to  a 
minimum  at  the  Brewster  angle  and  arg  R  changes  from  180*  to  0*.  Appendix  4 
gives  some  plots  of  the  reflection  coefficient  for  various  types  of  reflecting 
surface . 

The  above  formulae  assume  "smooth"  reflecting  surfaces.  At  the  wavelengths 
under  consideration  most  surfaces  will  appear  smooth  at  small  angles  of 
incidence,  but  at  larger  angles  diffraction  effects  will  become  apparent.  These 
effects  are  not  included  in  the  propagation  model . 


2.2  Calculation  of  path  loss. 

The  path  loss  between  two  antennas  of  gains  gT  and  gR,  at  a  frequency 
f  and  a  distance  d  apart  is  given  by, 


L-10  log  gT  +  10  log  gR  +  10  log  (c2/4s)  -  20  log  (f)  +  10  log  (l/4wd2) 


The  last  term  is  the  free  space  loss,  the  remaining  terms  take  account  of 
the  antenna  gains  and  the  aperture  of  the  receiving  antenna  which,  of  course, 
is  frequency  dependent. 

Fig. 2.1  illustrates  the  geometry  of  a  typical  multipath  situation.  If  the 
contributions  from  the  surface  vave  and  any  remote  structures  are  ignored  then 
the  signal  arriving  at  the  receiving  antenna  is  primarily  composed  of  the 
direct  ray  and  the  ground  reflected  ray.  The  lengths  of  these  paths  are, 


2 


For  the  direct  ray,  D  -  J{ (hl-h2)2  +  d2 ) 
For  Che  reflected  rey,  R  -  J( (hl+h2)2  +  d2) 
and  Che  grazing  angle  Is, 


f 


arc Can  ( 


hl4h2 

d 


) 


TV  iwAst 

h - d - *1 


Fig. 2.1  Multipath  geometry. 


The  propagation  model  uses  these  relationships  to  work  out  the  relative 
amplitudes  and  phases  of  the  direct  and  reflected  rays  at  the  receiving 
antenna.  Remember  that  the  phase  angle  between  the  rays  Is  a  combination  of  the 
path  length  difference  and  the  phase  change  due  to  the  reflection  coefficient. 
The  phasor  sum  of  the  two  rays  is  plotted  as  site  attenuation  versus  distance 
where  site  attenuation  Is  defined  as , 


Site  attenuation 


Power  transferred  by  direct  connection 
Power  transferred  by  propagation  paths 


3-  Propagation  prediction  software. 

A  complete  listing  of  the  program  is  contained  In  Appendix  2.  It  Is 
written  In  HP  BASIC  5.1  and  the  COMPLEX  binary  must  be  resident.  The  user  Is 
prompted  to  enter  frequency,  range,  antenna  heights  and  gains,  polarisation  and 
the  ground  constants.  Vertical  scaling  limits  are  set  by  the  variables  Ymax  and 
Ymln  and  may  be  changed  to  suit  a  particular  simulation. 
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4. 


Measurement  of  propagation  lost 


4.1  inttafluEtlfla. 

The  aim  of  the  propagation  measurement  trial  was  primarily  to  validate  the 
results  produced  by  the  propagation  model.  However,  the  choice  of  antenna  types 
and  the  heights  and  distances  involved  was  to  a  large  extent  determined  by  the 
application. 

With  reference  to  the  definition  of  site  attenuation  given  above  it  will 
be  seen  that  a  measurement  of  this  quantity  will  require  knowledge  of  the 
losses  in  the  cables  connecting  the  transmitter  to  the  transmit  antenna  and  the 
receive  antenna  to  the  receiver.  This  would  be  accomplished  by  connecting  the 
cables  directly  together  at  the  antennas  and  measuring  the  amount  of  power 
transferred  at  the  frequency  of  Interest.  With  the  antennas  reconnected  and 
separated  by  a  known  distance  a  second  measurement  of  transmitted  power  is 
made.  The  site  attenuation  is  then  the  ratio  of  these  two  levels  or  in 
decibels, 


Site  attenuation  -  PI  -  P2 

where  PI  -  power  received  by  direct  connection  (dB) 
F2  -  power  received  via  antennas  (dB) 


In  practice,  the  measurement  of  PI  need  only  be  performed  once  at  each 
frequency  of  measurement  and  used  for  any  subsequent  antenna  separation  and/or 
antenna  type. 


4.2  Description  of  equipment  employed. 

The  transmit  antennas  were  in  a  fixed  position  mounted  on  a  tripod  (except 
for  the  log-periodic  antenna  which  was  clamped  onto  a  wooden  stake  driven  into 
the  ground) .  The  receive  antennas  were  mounted  on  a  trolley  which  was  moved  by 
a  motorised  pulley  arrangement  up  to  a  maximum  distance  of  50  metres.  A 
serrated  disc  attached  to  one  of  the  trolley  wheels  interrupted  a  light  beam 
which  caused  pulses  to  be  sent  along  a  cable  to  a  counter  circuit  inside  the 
equipment  cabin  where  the  measurement  software  converted  the  pulses  into 
distance  travelled  in  metres. 

For  convenience  a  large  proportion  of  the  measurement  procedure  was 
carried  out  under  computer  control.  A  Hewlett-Packard  9000/310  machine  was  used 
with  dual  3  1/2  inch  disc  drives  for  storage  of  programs  and  data.  The  RF 
measurements  were  made  using  an  HP8753A  Network  Analyser  which  has  a  maximum 
output  level  of  +2SdBm  and  a  lOOdB  dynamic  range.  Since  single  frequency 
measurements  were  being  made  the  analyser  was  set  to  zero  span  which  enabled 
readlnga  to  be  taken  fast  enough  for  the  required  spatial  resolution.  In  order 
to  ensure  a  good  signal  to  noise  ratio  a  1  Watt  amplifier  module  was  used  at 
the  output  port  of  the  network  analyser.  A  listing  of  the  software  used  to 
control  the  measurements  is  presented  In  Appendix  3. 

The  distances  between  the  cabin  housing  the  measuring  equipment  and  the 
antennas  necessitated  long  lengths  of  coaxial  feeder  cable.  The  loss  of  100m  of 
UR43  at  1GHz  is  about  70dB  which  would  make  meaningful  measurements  impossible. 
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Therefore  Hellax  cable  was  used  where  possible  and  UR67  where  insufficient 
lengths  of  Heliax  were  available . 

The  following  table  lists  the  various  combinations  of  antennas  used  in  the 
trial . 


TRANSMIT  AKIfiMA 

Chase  UHA  9105 
(dipole, bicone) 

Chase  UHA  9105 
(dipole) 

CLP  5130-2 
Log-periodic 


RECEIVE  ANTENNA 

Singer 

(dipole) 

1/4  wave  whip 

CLP  5130-2 
Log-periodic 


FREQ.  RANGE  (MHz) 
100,300,500,700,900 

400,600 

100,300,500,700,900 


The  procedure  for  taking  a  typical  measurement  would  be  as  follows; 

1)  Run  measurement  software. 

2)  Calibrate  network  analyser  at  the  required  frequency. 

3)  Measure  cable  losses  at  this  frequency. 

4)  Set  up  required  antennas. 

5)  Move  trolley  toward  transmit  antenna  and  stop  at  2m  distance. 

6)  Ensure  distance  counter  has  reset. 

7)  Start  software  taking  readings. 

8)  Start  trolley  moving  away  from  transmit  antenna. 

9)  Stop  trolley  at  50m 

10)  Store  data  on  disc. 


4-3  DtsctivtimM  .trial  slit. 

All  the  propagation  trials  were  carried  out  on  the  antennas  field  at  RSRE 
South  Site.  This  site  was  far  from  ideal  being  quite  cluttered  and  not 
particularly  flat.  However,  rather  than  use  a  more  remote  site  with  the 
attendant  problems  of  transport,  power  and  communication  it  was  decided  to 
tolerate  the  non-ideal  conditions.  Fig. 4.1  shows  a  plan  of  the  site. 

The  measuring  equipment  was  housed  in  a  shielded  cabin  with  mains  power 
taken  from  a  nearby  15A  power  point. 
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FIXEO  IRW'KITO'V 


5. 


Results  of  the  propagation  measurement  trial 


The  results  obtained  from  the  trial  are  presented  In  a  graphical  format  in 
Appendix  1.  Each  graph  shows  site  attenuation  In  dB  versus  distance  In  metres 
both  for  the  measured  data  and  also  the  output  of  the  prediction  software.  Ho 
measurement  of  the  ground  constants  was  undertaken  and  therefore  typical  values 
were  used  for  the  prediction.  Small  variations  In  the  constants  have  little 
effect  on  the  predicted  results  except  to  vary  the  sharpness  of  nulls. 

The  first  set  of  measurements  were  intended  to  provide  standard  results 
against  which  the  predictions  of  the  propagation  model  could  be  compared. 
Half-wave  dipoles  were  used  for  most  measurements  except  at  the  lowest 
frequency  (100  MHz)  where  a  biconlc  antenna  had  to  be  used  as  the  transmit 
antenna.  Fig.Al.l  to  Fig.Al.10  show  the  results  of  these  measurements.  In  the 
prediction  software  the  gains  of  the  antennas  were  specified  as  2.13dBi,  the 
gain  of  a  half-wave  dipole  relative  to  an  isotrope. 

The  next  set  of  measurements  (Flg.Al.il  to  Fig.A1.14)  show  the  measured 
and  predicted  path  losses  between  a  dipole  (1.7m)  and  a  1/4  wave  vertically 
polarised  whip  close  to  ground  level.  In  Flg.Al.il  the  dipole  was  vertically 
polarised  and  the  operating  frequency  was  400  MHz.  In  Fig.Al.12  the 
transmitting  dipole  was  horizontally  polarised  and  the  whip  vertically 
polarised  but  the  simulation  assumed  both  antennas  horizontally  polarised.  This 
plot  therefore  gives  some  idea  of  the  loss  due  to  cross  polarisation.  Fig.Al.13 
and  Fig.Al.14  repeat  these  results  but  at  600  MHz. 

The  final  measurements  were  performed  using  wide-band  log-periodic 
antennas.  Results  were  obtained  every  200  MHz  from  100  MHz  to  900  MHz.  The 
predicted  results  still  assume  dipole  gains  and  so  an  indication  of  the  gain  of 
the  log-periodic  antennas  can  be  obtained. 


6.  Discussion  of  trial  results. 

It  can  be  seen  by  referring  to  the  results  in  Appendix  1  that  the  measured 
and  predicted  results  generally  agree  quite  closely.  However,  there  are  some 
interesting  discrepancies,  most  notably  with  regard  to  the  position  of  the 
nulls.  This  is  thought  to  be  a  consequence  of  the  variation  in  height  of  the 
ground  along  which  the  trolley  was  moved,  small  changes  in  height  can  produce 
significant  changes  in  null  position.  This  error  would  also  account  for  some  of 
the  discrepancy  at  longer  ranges,  for  example  in  Fig.Al.l,  if  the  nulls  between 
10m  and  15m  coincided  then  the  error  approaching  50m  would  not  be  so  great.  The 
position  of  the  nulls  is  also  affected  by  the  angle  of  the  ground  reflection 
coefficient,  but  with  horizontal  polarisation  the  angle  is  almost  always  close 
to  180* .  It  is  also  possible  that  the  surface  wave  changes  null  position  but  at 
frequencies  greater  than  about  100  MHz  the  effect  will  be  very  small. 

It  is  interesting  to  note  that  vertical  polarisation  suffers  far  less  from 
deep  nulls  than  horizontal.  Appendix  4  shows  that  this  is  because  the  magnitude 
of  the  ground  reflection  coefficient  is  generally  less  with  vertical 
polarisation  than  that  with  horizontal. 

Discrepancies  in  the  far  field  could  in  part  be  due  to  variation  of  the 
dipole  gains  with  frequency.  The  predictions  from  the  model  assume  constant 
gains  of  2.13  dBl. 
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A  further  source  of  error  will  arise  from  the  rather  arbitrary  choice  of 
ground  electrical  parameters  for  the  predicted  results.  Investigation  of  the 
effects  of  changing  the  ground  parameters  in  the  model  has  revealed  that 
vertical  polarisation  at  low  frequencies  is  most  sensitive  in  this  respect. 
With  horizontal  polarisation,  the  most  noticeable  effect  is  on  the  depth  of  the 
nulls.  These  conclusions  are  in  agreement  with  a  similar  investigation  carried 
out  in  Ref . 1. 

Other  possible  sources  of  error  are,  reflections  from  nearby  buildings  and 
passing  vehicles,  signal  cable  leakage  and  additional  radiation  due  to  braid 
currents  (imperfect  balanced  to  unbalanced  transformation). 


7.  Conclusions. 

This  report  has  presented  the  results  of  a  comparison  between  measured 
propagation  loss  and  predicted  values  from  a  simple  computer  model.  Despite  the 
simplicity  of  the  model  it  has  been  shown  to  agree  quite  closely  with  the 
measured  results.  Furthermore,  predictions  from  the  model  also  tally  with 
results  from  other  models  (Ref.l  and  Ref. 2). 

It  is  intended  to  incorporate  the  model  into  other  more  complex 
simulations.  This  Is  relatively  straightforward  since  the  prediction  program  is 
quite  short.  Data  generated  by  the  package  can  be  stored  in  arrays  for  further 
processing  rather  than  being  displayed  directly. 
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APPENDIX  1 


Results  of  propagation  loss  measurements . 
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Frequency:  700  MHz  Horizontal  polarisation 

TX:  1.7  ■  RX:  1.87  a  Conductivity:  .01  Er:  10 

- Predicted  -  Measured 


FIC.Al.l 
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Range (a) 

FI lenaae:RUN10  11  Rug  1989  11:17:05 

700MHz  HORIZONTAL  1.87a  MOIST  GROUND  CHRSE/SINGER 


Frequency:  700  HHz  Vertical  polarisation 

TX:  1.7  a  RX:  1.87  m  Conductivity:  .01  Er:  10 

-  Predicted  -  Measured 


FIG.A1.2 


(BP)  uoi*tnuen«  »*IS 


u 


Range (m) 

Ft  I an we i RUN 12  U  Rug  1989  11:38:27 

9 00MHz  VERTICAL  1.87m  MOIST  GROUND  CHRSE/SINGER 


Frequency:  900  MHz  Horizontal  polarisation 
TX:  1.7  ai  RX:  1.97  ■  Conductivity:  .01  Er 
- Predicted  -  Measured 


Fl lenaae:RUN13  11  Rug  1999  11:44:02 

900MHz  HORIZONTAL  1.97*  MOIST  GROUND  CHRSE/SINGER 


Frequency t  SOB  MHz  Horizontal  polarisation 


FIC.A1.5 


Rang* (a) 

Fl l*naa*:RUN16  11  Rug  I9B9  14:16:19 

500MHz  HORIZONTAL  1.67a  MOIST  GROUND  CHRSE/SINGER 


Frequency :  580  MHz  Vertical  polarisation 

TX:  1.7  ■  RXs  1.87  ■  Conductivity:  .01  Er:  10 

-  Predicted  -  Measured 


FIG.A1.6 


VERTICAL  1.87*  MOIST  GROUND  CHRSE/SINGER 


Frequency:  300  MHz  Verticil  polerlsetfon 


FIC.A1.7 
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Renge(m) 

Fl  leneae:RUN17  11  Rug  1989  14:38:39 

300MHz  VERTICRL  1.87«  MOIST  GROUND  CHRSE/SINGER  DIPOLES 


Frequency:  300  MHz  Horizontal  polarisation 

TX:  1.7  a  RX:  1.87  a  Conductivity:  .01  Er:  10 

- Predicted  - Measured 


HORIZONTAL  1.87*  MOIST  GROUND  CHRSE/SINGER  DIPOLES 


FIG.Al.9 


Frequency:  10B  MHz  Horizonte!  polerteetton 

TX:  1,7  c  RX:  1.87  •  Conductivity:  .01  Er:  10 

-  Predicted  -  Measured 


FIG.A1.10 


Frequency:  480  MHz  Verticil  potertsation 

TX:  1.7  ■  RX:  .1  ■  Conductivity:  .01  Er:  10 

- -  Predicted  -  Measured 


FIG.A1.11 


1/4  WHIP  VERT ICR.  1.87m/. In  DflMP  GROUND.  RAIN  CHRSE  DIPOLE 


Frequency:  400  MHz  Horizonte)  po Ur  I set  I on 

TX:  1.7  a  RXt  .1  •  Conductivity:  .01  Er:  10 

- Predicted  - Measured 


I 


FIC.A1.12 


1 


Frequency:  6B0  MHz  Vertical  polarisation 

TX:  1,7  i  RX:  .1  *  Conductivity:  .01  Er:  10 

-  Predicted  — - - —  Measured 


FIC.A1.13 


(OP)  uoiitnuenv  e*is 
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Range (aO 

Ft 1enaae:RUN23  14  Rug  1989  11:38:38 

600MHz  1/4  HRVE  HHIP  VERTICRL  C0. 1m>  CHRSE  DIPOLE  VERTICRL  (1.7a.)  DRMP 


Frequency:  600  MHz  Horizontal  polarisation 

TX:  1.7  ■  RX:  .la  Conductivity:  .01  Er:  10 

-  Predicted  -  Measured 


FIG.A1.14 


(8P)  uouanuenv  *t»S 
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Range(a) 

Ft  I en see: RUNS 4  14  Rug  19B9  11:44:50 

600MHz  1/4  WHIP  VERTICHL  CHRSE  DIPOLE  HORIZONTAL  (1.7a) 


Range (m) 

FI  leni«*:RUN25  14  Hug  1989  14:45:51 

100  HHz  HORIZ  HEIGHT  1.35  m  LOG  PERIODICAL OG  PERIODIC  DHMP  GROUND 


Frequency:  300  MHz  Horizonte!  polarisation 

TX:  1.36  a  RX:  1.36  a  Conductivity:  .01  Er:  10 

- Predicted  -  Meeeured 


FIC.A1.16 
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Frequency:  500  MHz  Horizontal  polarisation 

TX:  1.36  •  RX:  1.36  ■  Conductivity:  .01  Er:  10 

- Predicted  -  Measured 


FIC.A1.17 
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Range (n) 

Ft lena«e:RUN27  14  Rug  1989  15:04:26 

500  MHZ  HORIZ  1.35  m  RRINY/ORMP  LP/LP 


Frequency:  700  MHz  Horizontal  polarisation 

TX:  1.36  ■  RX:  1.36  a  Conductivity:  .01  Er:  10 

- Predicted  -  Measured 


FIG.A1.18 
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IZONTRL  LOG  /  LOG  1.35* 


Frequency:  900  MHz  Horizontal  polarisation 


Range (») 

Fi 1enaae:RUN29  14  Rug  1989  15:47:33 

900  MHz  HORIZ  1.36  m  DRMP  LP/LP 


APPENDIX  2 


Propagation  prediction  software 


29 


10  . . . . . . . . . 

20  t 
30  ! 

40!  SHORT  RANGE  VHF/CHF  PROPAGATION 
50!  PREDICTION  PROGRAM 

60  ! 

70  ! 

80  !  P.R. Bellamy  May  1989 

90  I  VER  2  August  1989 

100! 

110!  This  software  requires  BASIC  5.1  vlth  COMPLEX  binary 
120! 

130! . ■■■.■■. . . . . . . . . •••••••■••• 

140  COMPLEX  N2 ,R,B,A 

150  REAL  Free_space , Fresne 1 , D , Resultant , Freq , Range , HI , H2 , Er , Sigma , Inc 

160  REAL  Ymax , Ymin , L , C , EO , Lambda , Antenns_size , Antenna_f actor , W , Topline , Nextlin 

e 

170  REAL  Graze, Theta, Angle, Ref  coeff ,P1,P2,X,Y 

180  GIN1T 

190  KEY  LABELS  OFF 

200  DEG 

210  CLEAR  SCREEN 

220  PRINT  "  SHORT  RANGE  PROPAGATION  MODEL 


230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 


INPUT  "FREQUENCY  (MHz)?”, Freq 
INPUT  "MAX.  RANGE  (m)?", Range 
INPUT  "TX  ANTENNA  HEIGHT  (m)?\Hl 
INPUT  "RX  ANTENNA  HEIGHT  (m)?",H2 
INPUT  "TX  ANTENNA  GAIN  (dBi)?",Gt 
INPUT  "RX  ANTENNA  GAIN  <dBi)?",Gr 
INPUT  "POLARISATION  (H , V, Both) " , Pol$ 

INPUT  "RELATIVE  PERMITTIVITY?" ,Er 

INPUT  "CONDUCTIVITY?" .Sigma 

LINPUT  "OUTPUT  TO  CRT  OR  PLOTTER? ", Dev$ 

! 

IF  Dev$[l]-"P"  OR  Dev$[l]-"p"  THEN 
LINPUT  "PLOTTER  ADDRESS?" ,Add$ 

PLOTTER  IS  VAL(Add$) , "HPGL" 

ELSE 

PLOTTER  IS  CRT, "INTERNAL" 

END  IF 

Inc-Range/300  !  300  HORIZONTAL  POINTS 

! 

! 

1 

Ymax-100 

Ymln-0 

VIEWPORT  10,120,7,85 
PEN  1 

WINDOW  0, Range, Ymax, Ymin 
GRID  Range/10, 10,0, Ymax 
CLIP  OFF 
CS1ZE  3,-5 

MOVE  Range/2, Ymax+.07*(Ymax-Ymin) 

LORG  5 

LABEL  "Range (m)" 


550  FOR  1^0  TO  Rang*  STEF  Range/10 
560  MOVE  L,Ymax+.03*(Yaax-Yain) 

570  LABEL  L 

580  NEXT  L 

590  MOVE  -.08*Range,Yaax-<Ynax-YBin)/2 

600  LDIR  90 

610  LABEL  "Site  attenuation  (dB)" 

620  LDIR  0 

630  FOR  >Yiin  TO  Yaax  STEP  10 

640  MOVE  - . 03*Range,L 

650  LABEL  L 

660  NEXT  L 

670  ! 

680  ! 

690  I 

700  ! 

710  ! 

720  C«3 . E+8  I  SPEED  OF  LIGHT 

730  EO-8.854E-12 

740  Laabda-300/Freq 

750  Antenna_aize-Laabda/21  ROUGH  SIZE  OF  AERIAL 

760  Fresnel-2*Antenna_aize*Antenna  size/Laabda!  NEAR  FIELD  /  FAR  FIELD 
770  W-2*PI*Freq*l . E+6 

780  Antenna_f actor-Gt+Cr+10*LGT ( C*C/ (4*PI ) ) - 20*LGT ( Freq*l . E+6) 

790  ! 

800  ! 

810  ! 

820  Topllne-Ynin- . l*(Ymax-Ymin) 

830  Nextline-Ymin- . 05*(Ymax-Ymin) 

840  LORG  2 

850  MOVE  0, Topline 

860  LABEL  "Frequency ;Freq; "  MHz" 

870  ! 

880  ! 

8P3  MOVE  Range/2 .Topline 

900  IF  Pol$-"H"  THEN  LABEL  "Horizontal" 

910  IF  Pol$-"V"  THEN  LABEL  "Vertical" 

920  IF  Pol$-"B"  THEN 

930  PEN  2 

940  MOVE  Range/2 .Topline 

950  DRAW  Range/1 . 8 . Topline 

960  PENUP 

970  PEN  1 

980  LABEL  "  Horizontal" 

990  MOVE  Range/1. 4, Topline 

1000  PEN  3 

1010  DRAW  Range/1 . 3 .Topline 

1020  PENUP 

1030  PEN  1 

1040  LABEL  "  Vertical" 

1050  END  IF 
1060  I 

1070  ! 

1080  MOVE  0, Next line 

1090  LABEL  "TX  :";Hl;"a  RX  :";H2;’a" 

1100  MOVE  Range/2, Nextline 

1110  LABEL  "Conductivity: ";Sigaa;"  Er:";Er 

1120  ! 
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1130  I 
1140  ! 


1150  Next_pol : IF  Pol$-"V"  THEN 


1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 


PEN  3 
ELSE 
PEN  2 
END  IF 

I 

FOR  I>-Fre*nel  TO  Rang*  STEP  Inc 

P1-SQR( (HI - H2 ) * (HI - H2 ) +D*D )  I  DIRECT  PATH  LENGTH 

P2-SQR((H1+H2)*(H1+H2)+D*D)  1  REFLECTED  PATH  LENGTH 

Craze-ATN((H1+H2)/D)  I  GRAZING  ANGLE  IN  DEGREES 

Theta-360*((P2-Pl)/Laitbda-INT(  (P2 -PI) /Lambda))  I  PATH  DIFFERENCE  IN 


DEGREES  (l.e.  PHASE) 

J 

1**1 

CALCULATE  REFLECTION  COEFFICIENT 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

N2-CMPLX(Er , -Sigma/(W+E0) ) 

B— SQR(N2 - COS (Graze) *C0S (Graze ) ) 

A-SIN(Graze) 

IF  Pol$-"V“  THEN 
R-(N2*A-B)/(N2*A+B) 

ELSE 

R-(A-B)/(B+A) 

END  IF 

Ref_coef f-ABS (R) 

Angle-ARG(R) 

!************************************************************************* 

! 

! 

!  ADD  PHASORS . 

Free_space_d-SQR(30/(P1*P1))!  FREE  SPACE  LOSS:  DIRECT  RAY 
Free_space_r-SQR(30/(P2*P2) ) !  FREE  SPACE  LOSS:  REFLECTED  RAY 
X-Free_space_d+Free_space_r*Ref_coeff*COS(Theta+Angle) 
Y-Free_space_r*Ref_coeff*SIN(Theta+Angle) 

Resultant-(X*X+Y*Y)/(120*PI) 


PLOT  THE  RESULT. . . 

PLOT  D, -10*LGT (Resultant) -Antenna  factor 


NEXT  D 

IF  Pol$-"B"  THEN 
Pol$-’V" 

PENUP 

GOTO  Next_pol 
END  IF 
I 
! 

PENUP 

I 

END 
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1710  ! 

1720  ! 

1730  ! 

1740  I 

1750  1 ***********************************************************************«»* 


1760  ! 

TYPICAL  GROUND  CONSTANTS 

1770  ! 

1780  1 

SIGMA  (Mhos/M) 

EPSILON 

1790  ! 

Sea  water 

5 

80 

1800  ! 

Freah  water 

8e-3 

80 

1810  ! 

Dry  sandy  coastal 

2e-3 

10 

1820  ! 

Marshy  forested 

8e-3 

12 

1830  1 

Rich  agricultural 

le-2 

15 

1840  ! 

Pastoral  rolling 

5e-3 

13 

1850  ! 

Rocky,  nountainous 

2e-  3 

5 

1860  1 

Mountainous 

le-3 

5 

1870  ! 

Cities,  residential 

2e-3 

5 

1880  1 

Cities,  industrial 

le-4 

3 

1890  ! 
1900  I 
1910  I 
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Trial  control  software 


e 

‘j 

i 


1 


t 


J 


10 
20 
30 
40 

50 

51 

52 

53 
70 
80 

81  Scare:  t 

OPTION  BASE 
DEG 
I 
! 


SHORT  RANGE  PROPAGATION 


CONTROL  SOFTWARE 


TRIAL 


(c)  P.R. Bellamy  /  I. R. Dixon  1  Aug  1989 


90 

100 

110 

120 

130 

140 

150 

160 


Variable  declarationa 


! 

J  i****^*************^************^*^^**********^********************* 
!  MAIN  PROGRAM 

!*★*■**■*★*★***★****★★★*★★★★★***★★*★*★*■*•*******★*★*★★★*****★*★★★*★*★***★★*** 


DIM  Message${40] ,Answer$[l] ,File${10] ,K$[1] 

REAL  Prop_data(32767 ,2) 

INTEGER  Big_cog_teeth, Ht_cog_teeth ,Hax_readlng* ,Max_distance , Pulae_per_re 
v,Net_an,Gpio .Readings 

170  INTEGER  Cal_done .True .False ,Wanted_freq 

180  REAL  Freq , Dist_per  pulse, Gear_ratio,Wheel_circua,Max_dist_pulse. Null 
190  ! 

200 
210 
220 

230 

231 

240 

241 

250 

251 
260 
262 

270 

271 
280 
281 
282 

284 

285 

286 

287 

288 

290 

291 

300 

301 

303 

304 

306 

307 

308 

310 

311 

312 


GOSUB  Initial_yalues 
1 

lnput_freq(Freq) 

i 

S e t_up_8 7 5 3 ( Fr eq ) 

! 

Calibrate_8753 

! 

Check_freq( Freq, Wanted  freq) 

IF  NOT  Wanted_freq  THEN 
Entry_beep 

DISP  "CONTINUE  WITH  ";Freq;"  MHz  (y)?" 
Wait_key(K$) 

IF  K$0"Y"  AND  K$0"yn  THEN  GOTO  Start 
END  IF 
! 

Reset  counter(Gpio) 

! 

GOSUB  Take_readings 
1 

I GOSUB  Test_data 
! 

KEY  LABELS  OFF 

REDIM  Prop  data(Readlngs-l,2) 

I 

Plot_deta(Prop  data(*),Dist  per  pulse) 

I 

DISP  "Store  data  on  disc?  (n)" 
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313 

315 

316 

317 

318 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

333 

334 

335 
340 
350 
360 
370 
380 
I  0 
390 


Entry_beep 
Wait  key(K$) 

IF  K$-"N"  OR  K$-"n"  THEN 

DISP  "NO!  ARE  YOU  ABSOLUTELY  SURE?" 
Entry  beep 
Walt  key(K$) 

IF  K$-"Y"  OR  K$-"y"  THEN 
CLEAR  SCREEN 

PRINT  "PROGRAM  TERMINATED" 

PAUSE 
END  IF 
END  IF 
! 


Store_data(Prop  data(*)) 

I 

CLEAR  SCREEN 
DISP  "FINISHED." 

I************************************************************************* 


STOP 

Initial_values:  I 

PRINTER  IS  CRT 
CLEAR  SCREEN 
KEY  LABELS  OFF 

PRINT  CHR$(140) ; CHR$(132) ; *  SHORT  RANGE  PROPAGAT 

N  TRIAL  " ;CHR$(136) ;CHR$(128) 

PRINT  " 


Initialising  variables . 

400  GINIT 

410  Max_readings-5000  !Max  number  of  readings  (<32768) 

420  Max_distance-100  ! Distance  in  metres 

430  ! 

440  Wheel_circum-. 776 
450  Big_cog_teeth-162 
460  Lit_cog_teeth-20 

470  Gear_rat io-Big_cog_tee th/Li t_cog_teeth 

480  Pulse_per_rev-100 

490  Dist_per_pulse-Wheel_circura/(Gear_ratio*Pulse_per_rev) 

500  Max_dist  pulse-INT(Max_distance/Dist_per_pulse) -  32768 
510  ! 

520  Net  an-716 
530  Gpio-12 
540  ! 

550  White-1 
560  Red-2 
570  Yellow-3 
580  Green-4 
590  Cyan-5 
600  Blue-6 
610  Magenta-7 
620  ! 

630  PRINT  "done." 

640  RETURN 
650  ! 

660  ! 

670  Takejreadings :  ! 

680  CLEAR  SCREEN 
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681 

PRINT  " 

DATA  GATHE 

682 

KEY  LABELS  ON 

683 

CONTROL  2,2;1 

684 

ON  KEY  1  LABEL  " 

ABORT"  GOTO  Abort 

686 

PRINT  " 

PRESS  KEY  1  TO  ABORT 

MEASUREMENT” 

687 

! 

688 

Readings-1 

689 

Entry  beep 

691 

DISP  "PRESS  RETURN  TO  START 

MEASUREMENT. . 

692 

Walt  key(K$) 

693 

DISP  "TAKING  DATA . 

694 

REPEAT 

695 

OUTPUT  Net  an; "OUTPMARK; " 

696 

ENTER  Net_an;Prop_data (Readings, 2) .Null 

697 

CONTROL  Gpio,2;0 

698 

STATUS  Gpio , 3 ; Prop_data (Readings , 1 ) 

699 

CONTROL  Cpio,2;l 

700 

Readlngs-Readings+1 

701 

UNTIL  Readings>Max  readings 

702 

Err(3) 

703 

RETURN 

704 

Abort :  ! 

705 

Err(2) 

706 

RETURN 

920 

< 

921 

1 

922 

Test  data:  ! 

930 

CLEAR  SCREEN 

931 

PRINT  "  DATA 

GATHE! 

932 

PRINT  " 

ROUTINE" 


ROUTINE” 


WARNING:  THIS  ROUTINE  DOES  NOT  READ  THE  DISTANCE  COUNTER" 

934  KEY  LABELS  ON 

935  CONTROL  2,2 ;1 

936  ON  KEY  1  LABEL  "  ABORT"  GOTO  Abort2 

937  PRINT  » 


PRESS  KEY  1  TO  ABORT  MEASUREMENT" 

938  I 

939  Readlng8-1 

940  Entry_beep 

942  DISP  "PRESS  RETURN  TO  START  MEASUREMENT..." 

943  Walt_key(K$) 

944  DISP  "TAKING  DATA . 

945  TO-TIMEDATE 

946  REPEAT 

947  OUTPUT  Ne  t_an ; "OUTPMARK ; " 

948  ENTER  Net_an;Prop_data(Readings , 2) .Null .Null 

949  Prop_data (Readings , 1)-Readings*30 . 0 

950  Readlngs-Readings+1 

951  UNTIL  Readings>Max  readings 

952  Err(3) 
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953 

RETURN 

954 

Abort2 :  1 

955 

Tl-TIMEDATE 

956 

PRINT  "TIME  PER  READING:  (T1 -TO) /Readings ;  "  SECS" 

957 

Err(2) 

958 

RETURN 

959 

! 

960 

! 

961 

! 

962 

END 

963 

964 

!  SUBPROGRAM 

DEFINITIONS 

965 

Set  up  8753:  I 

966 

SUB  Set  up  8753 (REAL  Freq) 

967 

Na-716 

968 

ABORT  7 

969 

OUTPUT  Na; "PRES ; " 

I  PRESET  8753 

970 

OUTPUT  Na; "FORM3 ; POIN  3;" 

! FORMAT  3,  3  POINTS 

971 

OUTPUT  Na; "LOGM; " 

! LOG  MAG  DISPLAY 

972 

OUTPUT  Na; "CHAN1; " 

! CHANNEL  1  ACTIVE 

975 

OUTPUT  Na; "IFBW30 ; " 

1 1F  BANDWIDTH 

976 

OUTPUT  Na; "POWE  0;" 

! SOURCE  POWER 

977 

OUTPUT  Na; "S21 ; " 

(TRANSMISSION  MEASUREMENT 

978 

OUTPUT  Na; "CENT  ";Freq;"  MHZ; SPAN 

0  MHZ;"  (SET  CENTRE  FREQ  AND 

979 

! 

980 

A$— CHR$  (34)6."  PROPCAL"  &CHR$  ( 34 ) 

(TITLE  REG  5 

981 

OUTPUT  Na ;  "TITR5 ;  "6>A$6i"  ;  " 

982 

SUBEND 

983 

1 

984 

! 

985 

Calibrate  8753:  ! 

986 

SUB  Calibrate  8753 

987 

CLEAR  SCREEN 

988 

PRINT  "  8  7 

53  CALIBRATION 

989 

PRINT  " 

990 

Entry  beep 

991 

INPUT  "LOAD  CAL  DATA  FROM  8753  REG 

5?  (Y/N)"t  Reply? 

992 

IF  Reply$-"Y"  OR  Reply$-"y"  THEN 

993 

OUTPUT  716; "OPC? ;RECA5; " 

994 

ENTER  716;Reply 

995 

ELSE 

996 

OUTPUT  716 ; "CALKN50; " 

997 

OUTPUT  716 ; "CALIFUL2 ; " 

(FULL  2  PORT  CALIBRATION 

998 

OUTPUT  716 ; "REFL; " 

999 

1 

1000 

DISP  "CONNECT  OPEN  TO  PORT  1  THEN 

PRESS  RETURN. . 

1001 

CALL  Walt  key(K$) 

1002 

OUTPUT  716; "C1ASS11A;0PC? ;STANB; " 

1003 

ENTER  716; Reply 

1004 

! 

1005 

DISP  "CONNECT  SHORT  TO  PORT  1  THEN  PRESS  RETURN. . 
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1006  CALL  Walt  key(K$) 

1007  OUTPUT  716; "CLASS11B;0PC? ;STANB; " 

1008  ENTER  716; Reply 

1009  I 

1010  DISP  "CONNECT  LOAD  TO  PORT  1  THEN  PRESS  RETURN. 

1011  CALL  Wait  key(K$) 

1012  OUTPUT  716 ; "CLASS11C;0PC? ;STANA; " 

1013  ENTER  716 ;Reply 

1014  ! 

1015  DISP  "CONNECT  OPEN  TO  PORT  2  THEN  PRESS  RETURN. 

1016  CALL  Wait  key(K$) 

1017  OUTPUT  716; "CLASS22A;OPC? ;STANB; " 

1018  ENTER  716;R#ply 

1019  ! 

1020  DISP  "CONNECT  SHORT  TO  PORT  2  THEN  PRESS  RETURN 

1021  CALL  Wait_key(K$) 

1022  OUTPUT  716; "CLASS22B;OPC?;STANB;" 

1023  ENTER  716;Reply 

1024  ! 

1025  DISP  "CONNECT  LOAD  TO  PORT  2  THEN  PRESS  RETURN. 

1026  CALL  Wait_key(K$) 

1027  OUTPUT  716; "CLASS22C;OPC? ; STANA; " 

1028  ENTER  716;Reply 

1029  OUTPUT  716; "REFD; " 

1030  I 

1031  I 

1032  DISP  "CONNECT  THRU  THEN  PRESS  RETURN..." 

1033  CALL  Wait  key(K$) 

1034  OUTPUT  716; "TRAN;" 

1035  ! 

1036  OUTPUT  716 ; "OPC? ; FVDT ; " 

1037  ENTER  716;Reply 

1038  ! 

1039  OUTPUT  716 ; "OPC? ;REVT; " 

1040  ENTER  716; Reply 

1041  ! 

1042  OUTPUT  716 ; "OPC7 ; FWDM ; " 

1043  ENTER  716; Reply 

1044  ! 

1045  OUTPUT  716; "OPC? ;REVM; " 

1046  ENTER  716;Reply 

1047  OUTPUT  716; "TRAD;" 

1048  ! 

1049  OUTPUT  716; "ISOL; " 

1050  OUTPUT  716;"0MII;" 

1051  OUTPUT  716;"ISOD;" 

1052  ! 

1053  OUTPUT  716;"SAV2;" 

1054  OUTPUT  716 ; "OPC? ; SAVE5 ; ” 

1055  ENTER  716; Reply 

1056  Cal_done-True 

1057  END  IF 

1058  DISP  "CALIBRATION  COMPLETE" 

1059  WAIT  1 

1060  CLEAR  SCREEN 

1061  SUBEND 

1062  ! 

1063  ! 
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1064  Wait_key:  f 

1065  SUB~Wa i t_key (K$ ) 

1066  DIM  T$[10] 

1067  ON  KBD  GOTO  Quit 

1068  LOOP 

1069  END  LOOP 

1070  Quit: ! 

1071  T$-KBD$ 

1072  K$-T$ [1,1] 

1073  SUBEND 

1074  ! 

1075  1 


1076  Input_freq:  I 

1077  SUB  Input_f req (REAL  Freq) 

1078  Entry_beep 

1079  LOOP 

1080  INPUT  "Frequency  for  this  measurement  (MHz)  ",Freq 

1081  EXIT  IF  NOT  (Freq<100)  AND  NOT  (Freq>1000) 

1082  Err(l) 

1083  END  LOOP 

1084  SUBEND 

1085  ! 

1086  ! 

1087  Err:  !  MASTER  ERROR  HANDLING  ROUTINE 


1088 

1089 

1090 

1091 

1092 

1093 

1094 

1095 

1096 

1097 

1098 

1099 

1100 
1101 
1102 

1103 

1104 
1106 

1107 

1108 

1109 

1110 
1111 
1112 

1113 

1114 

1115 


SUB  Err (INTEGER  Err_num) 

SELECT  Err_num 

j 

CASE  1 

DISP  "FREQUENCY  OUT  OF  RANGE" ; 

! 

CASE  2 

DISP  "MEASUREMENT  ABORTED"; 

! 

CASE  3 

DISP  "MAXIMUM  NUMBER  OF  READINGS  EXCEEDED" ; 


CASE  4 

DISP  "CALIBRATION  DATA  NOT  VALID  AT  THIS  FREQUENCY" 

! 

CASE  5 

DISP  "RESET  FAILED” 

i 

CASE  ELSE 
END  SELECT 

DISP  "  <SPACE>  TO  CONT" 

BEEP  1000,. 2 
Walt_key(K$) 

DISP  ■ 

SUBEND 

! 

! 


1116  Entry_beep :  ! 

1117  SUB  Entry_beep 

1118  BEEP  900 ,.05 

1119  SUBEND 

1120  I 

1121  ! 


1122  Reset  counter:  ! 
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1123  SUB  Reset_counter(INTEGER  Gpio) 

1124  INTEGER  A 

1125  1  Output*  a  reset  pulse  to  the  distance  counter  on  CTL1 

1126  Do_agaln:  ! 

1128  CONTROL  Gpio, 2 ;1 I  CTL1-1  (RESET  HIGH),  CTLO-O  (CONTROL  LOW) 

1129  WAIT  .1 

1130  CONTROL  Gpio,2;3!  CTL1-0  (RESET  LOW  ),  CTLO-O  (CONTROL  LOW) 

1131  WAIT  . 1 

1132  CONTROL  Gpio,2;l!  CTL1-1  (RESET  HIGH),  CTLO-O  (CONTROL  LOW) 

1133  I 

1134  DISP  "RESET  PULSE  SENT  -  MOVE  WHEEL  THEN  PRESS  A  KEY. 

1135  Walt  key(A$) 

1137  CONTROL  Cplo,2;0 

1138  STATUS  Gpio,3;A 

1139  CONTROL  Gplo,2;l 

1140  IF  A<0  OR  A>100  THEN 

1141  Err(5) 

1142  GOTO  Do  again 

1143  END  IF 

1144  SUBEND 

1145  ! 

1146  1 

1147  Check  freq:  ! 

1148  SUB  Check_f req (REAL  Freq, INTEGER  Flag) 

1149  REAL  F 

1150  Flag-1 

1151  OUTPUT  716; "CENT?;" 

1152  ENTER  716 ;F 

1153  IF  FOFreq*l .  E+6  THEN 

1154  CALL  Err(4) 

1155  Freq-F/1 . E+6 

1156  Flag-0 

1157  END  IF 

1158  SUBEND 

1159  ! 

1160  ! 

1161  Plot_data:  ! 

1162  SUB  Plot_data(REAL  Prop_data(*) ,Dist_per_pulse) 

1163  INTEGER  White.Red, Yellow, Green, Cyan, Blue, Magenta 

1164  INTEGER  I 

1165  Whlte-1 

1166  Red-2 

1167  Yellow-3 

1168  Green-4 

1169  Cyan-5 

1170  Blue-6 

1171  Magenta-7 

1172  ! 

1173  CLEAR  SCREEN 

1174  GINIT 

1175  GCLEAR 

1176  GRAPHICS  ON 

1177  INPUT  "Start  distance: ".Start  pos 

1178  INPUT  "Aap  gain  (dB):",Gain 

1179  PEN  White 

1180  Readings-S IZE ( Prop  data.l) 

1181  I 

1182  VIEWPORT  10,120,28,95 
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1183  WINDOW  0,50,100,0 

1184  AREA  PEN  Blue 

1185  MOVE  0,0 

1186  RECTANGLE  50, 100, FILL 

1187  AXES  1,2,0,100,10,5,3 

1188  ! 

1189  CLIP  OFF 

1190  CSIZE  2.5 

1191  LD1R  0 

1192  LORG  8 

1193  FOR  1-0  TO  100  STEP  10 

1194  MOVE  0,1 

1195  LABEL  1 

1196  NEXT  I 

1197  CSIZE  3 

1198  LORG  4 

1199  LDIR  90 

1200  MOVE  -3,50 

1201  LABEL  "Sice  attenuation  (dB)” 

1202  ! 

1203  CSIZE  2 . 5 

1204  LDIR  0 

1205  LORG  6 

1206  FOR  1-0  TO  50  STEP  10 

1207  MOVE  I , 100 

1208  LABEL  I 

1209  NEXT  I 

1210  CSIZE  3 

1211  MOVE  25,102 

1212  LABEL  -Range  (m)" 

1213  ! 

1214  MOVE  40,0 

1215  LDIR  0 

1216  LORG  4 

1217  CSIZE  2.7 

1218  LABEL  DATES (TIMEDATE) ; "  TIMES (TIMEDATE) 

1219  ! 

1220  FOR  1-2  TO  Readings 

1221  Dist-FNDistance(Prop_data(I,l) ,Dist_per_pulse)+Start_pos 

1222  Prop  data(I,2)--Prop_data(I,2)+Cain 

1223  PLOT  Dist,Prop_data(I,2) 

1224  Prop_data( 1 , 1 ) -Dis  t 

1225  NEXT  I 

1226  ! 

1227  SUBEND 

1228  I 

1229  t 

1230  DEF  FNDi stance (REAL  Counter ,Dist_per_pulse) 

1231  IF  CountercO  THEN 

1232  RETURN  (Counter+65536)*Dist_per_pulse 

1233  ELSE 

1234  RETURN  Counter*Dist  per_pulse 

1235  END  IF 

1236  FNEND 

1237  1 

1238  ! 

1239  Store_data:  I 

1240  t  Creates  file  on  which  the  data  In  Prop_data  Is  stored 
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1241  !  Summary?  i*  a  160  character  (2  line)  summary 

1242  SUB  Store_data(REAL  Prop_data(*)) 

1243  ON  ERROR  GOTO  Trap 

1244  CLEAR  SCREEN 

1245  DIM  Sumaary$[160] 

1246  DIM  File?[25] 

1247  Readings-SIZE(Prop  data.l) 

1248  ! 

1249  ! 

1250  PRINT  TABXY(1,21); "SUMMARY:  (2  lines  max)" 

1251  LINPUT  Summary? 

1252  Inp : LINPUT  "File:mass  storage  unit?", File? 

1253  DISP  "Creating  file..." 

1254  CREATE  BOAT  File?,Readings*16+500,l  I  NOTE:  BASIC  5  seems  to  w 

ork  out 

1255  1  the  correct  file  size  itself! 

1256  ASSIGN  @Path  TO  File? 

1257  DISP  "Storing  data..." 

1258  OUTPUT  @Path ; Summary? .TIMEDATE .Readings , Prop_data(*) 

1259  ASSIGN  @Path  TO  * 

1260  DISP  " 

1261  SUBEXIT 

1262  ! 

1263  Trap:  ! 

1264  DISP  ERRM? ; "  <SPACE>  TO  CONT" 

1265  BEEP  1000,. 2 

1266  Wait_key(K$) 

1267  GOTO  Inp 

1268  SUBEND 
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APPENDIX  4 


Ground  reflection  coefficient  plots 
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Rrg(p)  ^  Mag(p)  i  i  Arg(p)  Hag(p) 
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Grazing  angle 


Relative  permittivity:  88  Conductivity:  5  mhos-'m 


Horizontal  polarisation 


L _ _ l _ ...  J _ I  - I  I - L  . 1 


Frequency:  l.E+8  Hz 
Relative  permittivity:  88 
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SEA  WATER 

Conductivity:  5  mhos/m 


Reletlve  permtttt vlty:  12  Conductivity:  .006  mhos/m 


Frequency:  l.E+8  Hz 
Reletlve  permittivity:  18 
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MRR5HY  FORESTED 
Conductivity:  .006  mhos/m 


(<*)•)««  0*)8jh  1  '  <B)B«w  (d)Bjb 


Frequency:  l.E+B  Hz  URBAN  RESIDENTIAL 

Relative  permittivity:  5  Conductivity:  .BBS  mhos/m 
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Abstract  This  report  describes  a  VHF/UHF  propagation  prediction  model  and  a 
series  of  field  measurements  aimed  at  verifying  the  results  of  the  model.. 
Comparisons  between  the  model  and  the  trial  results  were  made  at  spot  frequencies 
from  100  MHz  to  900  MHz  for  several  types  of  antenna  and  for  various  heights 
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The  predicted  path  losses  from  the  model  were  found  to  agree  quite  well  with  the 
results  of  the  field  measurements.  It  should  be  noted  that  the  prediction 
model  was  never  intended  to  be  a  comprehensive  tool,  merely  a  utility  package 
which  could  be  incorporated  into  other  software.  — v  fX_r  t  *  c,'"  -  .. 

Continued  of  next  page 


ssoms 


CONTINUED 


The  program  was  written  on  an  HP  9000  series  model  310  computer  using  BASIC 
5.1.  This  language  has  a  set  of  commands  which  allow  direct  manipulation  of 
complex  values,  consequently,  translation  into  another  language  or  version  of 


BASIC  not  having  these  facilities  will  require  extra  programming. 
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